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ABSTRACT 


NOAA-7 AVHRR data from April 1982 are used to demonstrate a method for 
building an aerosol climatology for the North and South Pacific Ocean. Sunglint 1s 
removed by deleting an area determined from visual inspection of many images. 
Red-visible (0.63um), near-infrared (0.864m) and infrared (1ljtm) channels, and the 
ratio of visible albedo to near-infrared albedo are used to eliminate areas of cloudiness. 
The ratio values are also used to indicate whether small (continental) or large (marine) 
aerosol particles predominate. Once a “clean” image of radiances due primarily to 
scattering by aerosols is obtained, satellite pixel values are averaged over one degree 
boxes from 50°N to 50°S, 1 10°E to 70°W. The mean visible and near-infrared albedos, 
mean ratio and standard deviations are computed and displayed as contoured fields 
over the North and South Pacific Ocean. These charts and individual examples are 


used to estimate spatial and temporal changes in aerosol concentration. 
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I. INTRODUCTION 


Atmospheric aerosols scatter electromagnetic radiation over a broad spectrum of 
wavelengths, including those ranging from the visible (0.4j1m) to infrared (121m). 
Aerosol particles range ın size from 0.00lum to 100um, however a much smaller range, 
0.1um to 108850. most effectively contributes to the scattering of the visible radiation. 
Many of the physical characteristics of scattered radiation are parameterized in terms 


of the size parameter, s: 


)1.1( << و 


For the wavelengths and particle sizes cited above, scattering by aerosols 1s significant 
for s values from approximately 0.05 to 160. The smallest s values occur with small 
particles and long wavelengths (Rayleigh scattering), whereas large s values (geometric 
optics for s> ~ 20) are obtained with large particles and short wavelengths. However, 
since most scattering occurs between the two extremes of particle size and wavelength, 
where s 1s on the order of one, analysis of the scattering requires application of Mie 
theory. 

Atmospheric aerosols, depending on their chemistry, shape and size distribution, 
can be divided into two primary groups. Continental particles are composed mainly of 
silicates and carbon; marine particles are composed mostly of sea salt. Although the 
distribution of particles in the atmosphere at any place and time may be a combination 
of these types, continental distributions typically have many more small particles 
(r< lium) than marine distributions. 

This thesis shows that an aerosol climatology can be built from satellite radiance 
data. Specifically, NOAA-7 Advanced Very High Resolution (AVHRR) data from 
April 1982 have been used to construct contoured fields of red-visible and near-infrared 
albedos, and the ratio of these albedos, as well as their standard deviations and the 
quantity of data used in the calculations. The procedures used to eliminate sunglint, 


land and rlouds are also described. 


A. BACKGROUND 

As early as 1970 it was recognized that dust from Saharan storms could be 
detected over the North Atlantic Ocean in satellite images (Shenk and Curran, 1974). 
These dust clouds were mapped in 1979 with isopleths of optical depth from NOAA-3 
data (Carlson, 1979). Less obvious background atmospheric aerosols were first seen as 
gray shades in Very High Resolution (VHR) visible images over water from Defense 
Meteorological Satellite Program (DMSP) satellites (Fett and Isaacs, 1979). Also, 
Earth Resources Technology Satellite One (ERTS-1) visible data were used to establish 
a linear upwelling radiance-aerosol optical depth relationship over water at 0.55pm, 
0.65pm, and 0.754m (Griggs, 1975). Optical depth was thus established to be one 
measure of aerosol concentration. Subsequent studies of data from different satellite 
sensors with higher spatial resolution have also confirmed that upwelling radiance and 
optical depth are nearly linearly related (Carlson and Wendling, 1977; Griggs, 1979; 
Durkee et al., 1986). Griggs (1983) concluded that the linear relationship can be used 
globally to monitor aerosol optical depth from satellite detected radiances. He also 
suggested that the aerosol size distribution could be inferred from a comparison of 
radiances from channel | (0.63um) and channel 2 (0.86um) of the AVHRR on 
NOAA-6. Takashima and Takayama (1983) used visible and infrared AVHRR data 
from NOAA-6 to estimate radiance variations resulting from sunglint and from ocean 
wave heights resulting from changing wind speed. In addition, they calculated 
radiation transmitted through the atmosphere in the visible and infrared spectral 
region, thereby accounting for absorption by atmospheric gases and scattering by 
aerosols. 

Durkee (1984) showed how aerosol particle size variations affect satellite detected 
upwelling radiance. Since marine aerosol particles are produced by evaporating 
droplets from bursting seawater bubbles (Woodcock, 1953) and are therefore composed 
primarily of water soluble material, particle size is significantly affected by relative 
humidity changes (Hânel, 1976). As particle size increases, extinction increases, which 
implies that upwelling radiance should increase. Durkee (1984) verified this reasoning 
for a case study in southern California. 

Meszaros and Vissy (1974) measured the concentration, size distribution and 
chemistry of aerosols over the Southern Ocean. Concentrations were less than those in 
the northern hemisphere for all sizes. The size distribution ranged between 0.03um and 
64.0um, with a maximum near 0.lftm. The particles were mainly sea salt and sulphur 


compounds. These results defined a background aerosol, undisturbed by land sources. 
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B. OBJECTIVES/MOTIVATION 

In this thesis the relationship between satellite detected upwelling radiance and 
optical depth is used to establish the physical basis for an aerosol climatology of the 
North and South Pacific Ocean. Brightness counts from channel I and channel 2 of 
the AVHRR are converted to albedos using calibration algorithms provided by NOAA 
(Kidwell, 1984). Combinations of these values are averaged in 1° squares over the 
entire basin. These average values are then used to build a data base for the North 
and South Pacific Ocean. The data used to demonstrate the method were AVHRR 
measurements from April 1982. 

A bias toward fair weather radiance values is anticipated because cloudy portions 
of the satellite pass have been excluded. Aerosol particles in cloudy regions, which 
have water-saturated air, act as condensation nuclei and thus the aerosol can grow 
rapidly to cloud droplet size. The size distribution of this study’s data sample may be 
different from that of the North and South Pacific Ocean, cloudy and cloud-free. The 
data must be interpreted accordingly. 

The ratio of albedo at channel 1 to albedo at channel 2 ıs also averaged in 1° 
squares over the North and South Pacific Ocean. Ratios of albedos at different 
wavelengths vary as the aerosol particle concentration varies. Durkee (1984) and 
Bulfinch (1986) showed that as the relative concentration of smaller particles increase, 
the ratio values of AVHRR channel 1 albedo/channel 2 albedo increase. Therefore 
horizontal variations in ratio values indicate horizontal variations in size distribution, 
with large ratio values indicating more small (r < lum) particles and a more continental 
distribution type. 

If the climatological background aerosol concentration over the North and South 
Pacific Ocean is known, it should be straightforward to detect changes in aerosol size 
distribution and composition resulting from an intrusion of continental aerosol 
particles into the troposphere over the Pacific Ocean. Iwasaka et al. (1983) studied an 
Asian dust storm which occurred in April 1979. Ground based measurements of the 
aerosol particles yielded a representative size distribution function and extinction 
coefficient. These data, as well as a Japanese Geosynchronous Meteorological Satellite 
visible image clearly showing the dust cloud, indicate that Asian dust events moving 
over the North Pacific Ocean are similar to Saharan dust storms which move over the 
North Atlantic Ocean (Carlson and Wendling, 1977; Carlson, 1979; Norton, et al., 
1980; Shenk and Curran, 1974). A case study will be presented which documents an 


p» 


event in the eastern North Pacific Ocean that may be related to the eruption of El 
Chichon, Mexico in March-April 1982. The ability to identify and track these clouds 
across the North Pacific Ocean may lead to improved forecasts of electro-optical 
system performance and improved visibility forecasts. The climatological data have 
important applications to climate studies due to the dust particles’ role in cloud physics 
and radiative transfer processes. Furthermore, this information bears directly on 
chemical oceanography because the dust can be carried over a large area and can 


contribute significantly to deep ocean sedimentation (Iwasaka, et al., 1983). 


C. ORGANIZATION 

Chapter 2 will discuss the theory of radiative scattering by atmospheric aerosols, 
the basic equations of radiative transfer and the assumptions applied to simplify them. 
The data set and acquisition equipment are described in Chapter 3, as 1s the procedure 
used to remove contaminants (sunglint and clouds). Since the objective 1s an aerosol 
climatology, it must be established that only radiation scattered by aerosols is being 
sensed by the instrument, not by other atmospheric constituents or phenomena. A 
case study of one area in the North Pacific Ocean is given, with an example of the 
contaminant discriminating process. Conclusions and recommendations for future 


study are stated in Chapter 4. 
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II. THEORY 


A. THE ELECTROMAGNETIC SPECTRUM 

For remote sensing of aerosols, we detect the interaction of visible and 
near-infrared radiation with aerosol particles dispersed in air. This is called the inverse 
problem; i.e., information about the particles doing the scattering 1s inferred from an 
analysis of the field of scattered radiation (Bohren and Huffman, 1983). 

The radiation detected by a satellite sensor has been scattered by air molecules, 
aerosol particles suspended in the air and the earth's surface. Molecular scattering is 
described by Rayleigh theory (s « 1) Molecular scattering makes up a significant 
portion of the total atmospheric scattering, but, as will be discussed below, does not 
vary significantly with position. At red visible and near-infrared wavelengths, 
reflectance by an ocean surface is very small (Ramsey, 1968). Therefore any variation 
in the radiation field is primarily due to aerosol scattering. Mie theory describes 
aerosol scattering (s ~ 1). The radiation field is affected by the size and number of 
particles present, and their chemistry. All of these vary significantly in the horizontal 


due to their sources. 


B. THE RADIATIVE TRANSFER EQUATION 

The radiative transfer equation (RTE) can be used to completely describe the 
radiation field at a satellite sensor,accounting for attenuation (absorption, scattering 
and reflection out of the field of view) and intensification (emission, and scattering and 
reflection into the field of view). 

1. Absorption and Scattering 

Atmospheric gases have absorption bands which, when combined, give the 

resultant absorption spectrum of the earth’s atmosphere (see Fig. 2.1). Some sensors, 
such as the human eye, are “designed” to take advantage of the “windows” where the 
solar radiance is absorbed least; that is, where the sun’s energy is able to penetrate the 
atmosphere. In the window regions the radiation can also be scattered back toward 
the sensor by gas molecules, cloud droplets, aerosols, or the earth’s surface (land or 
sea). AVHRR channel 1, at 0.55ptm to 0.68um, and channel 2, at 0.704m to 1.00pm, 


span two of the atmosphere’s absorption windows. 
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The complex refractive index represents both the scattering and absorbing 


ability of a substance. It 1s defined as: 
m=n+ iv. (2.1) 


The real part, n, is the ratio of the speed of light in one substance to that in another; 
the magnitude of scattering is partially governed by the magnitude of n. The imaginary 
part, v, is a measure of the absorptivity of a substance. Wavelength strongly affects v, 
but v is no larger than 1074 for wavelengths less than lfm (Shettle and Fenn, 1979); 
therefore, there is little absorption by aerosols of the radiation sensed by AVHRR 
channels 1 and 2. 

Scattering by aerosols is described using the scattering coefficient, p 


scat? 
extinction is given by the extinction coefficient, Bext (see Bohren and Huffman, 1983): 


Bscat = THO ٥٢٦ (2.2) 
and 

Bext = f tr? Qexe(mar){AN(0)/ar} dr, (2.3) 
where 


r = particle radius, 

m = complex refractive index, 
Q.cat = scattering efficiency, 

Qext = extinction efficiency and 
dN(r)/dr = particle size distribution. 


The scattering and extinction coefficients are related in the following way: 


)2.4( :سط 8 


Therefore the atmospheric scattering and extinction of aerosols is the vertical sum 
through the atmosphere of the product of three terms: the cross sectional area of the 
particles (rr), the scattering or extinction efficiency (Q.oa4, Rex)» and the size 


distribution (dN/dr). Fig. 2.2 illustrates the variation of these terms witn radius. The 


16 


relative magnitudes of these terms determines the magnitude of the extinction 


coefficient which in turn gives the optical depth, defined as: 


n foBextda (2.5) 


where H is the height of the satellite. It can be seen from Fig. 2.2 that when r is small, 
nr? and QM are small while dN/dr 1s large; when r is large, nr? Is large, and dN/dr 


and رم‎ become small. The result is that D... is affected by a bounded region of 


ext 
particle sizes. For visible wavelengths and a marine aerosol size distribution model at 
80% relative humidity this range is from approximately 0.54m to 5.0um (Durkee, 
1984). Because the imaginary part of the complex index of refraction is so small at 


these wavelengths (v ^ 1074), p This means that marine aerosols are 


scat ~ Pext: 
efficient scatterers. The ratio: 


EO 5 Bscat/Pext: (2.6) 


called the single scattering albedo, is approximately I in this case. 
2. Angular Dependence 

The scattering phase function, p($2,$2.), describes the angular dependence of 
scattering. If we assume the particles are spherical (necessary to get an exact solution 
to Maxwell's equations for electromagnetic waves, and one which is essential to keep 
the problem computationally tractable), Mie theory can be used to calculate p($2,X2..). 
For a complete description, see Bohren and Huffman, 1983. Fig. 2.3 shows the phase 
function calculated from Mie theory at 80% relative humidity. Marine aerosols can be 
assumed spherical if the relative humidity is high, because they are hygroscopic and 
therefore tend toward spheres as they gain water on their surface. For the marine 
aerosol at 80% relative humidity, particles primarily scatter forward (Shettle and Fenn, 
1979). 

The albedo of the ocean surface is nearly zero in the red-visual to 
near-infrared region of the spectrum (Ramsey, 1968), but sunglint 1s observed because 
specular reflection occurs off wave facets oriented such that incoming solar radiation is 
reflected into the sensor’s view. The size of the sunglint area increases as the wind 


velocity increases (Cox and Munk, 1954). 
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Fig. 2.3 Single scattering phase function for the model of marine 
particles at 80% relative humidity after Shettle and Fenn (1979). 
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Exe ee TIE IS C NUS C ۹ 


Due to the angular dependence of particle scattering and surface reflectance, 
the RTE depends on sun-earth-satellite geometry. Fig. 2.4 identifies all of the angles 
involved. These are important because the amount of radiation received by the 
satellite sensor depends on the view direction of the sensor. The intensity detected 
varies with the observing zenith (Û) and azimuth (@) angles. If the upwelling radiance 
depends on reflection of solar radiation, solar zenith (0) and azimuth (@,) must also 
be accounted for. The angle between the incoming solar beam and the beam received 
at the sensor (@) is a function of these four angles. The angular direction of the 
upwelling beam described by (0,@) is called £2 and that described by (8,,,) is £25. 

3. The General RTE and Simplifying Assumptions 

We can now write the RTE. A basic assumption is that the atmosphere is 
plane parallel; i.e., the top and bottom of the atmosphere are two parallel planes. If 
optical properties of the atmosphere are assumed to be horizontally homogeneous, the 


general form of the RTE is given by (see for example Liou, 1980): 


udL(t,02)/dt = L(7,02) (2.7) 
- (85/41) f gp L(7,2’)p(Q,Q’)dQ’ 
- (0,/40)TF pP(Q,-Q,)e Vo, 
where 
L = diffuse radiance, 


t — optical depth, 


©, — single scattering albedo, 
u = così, 
0۳٣ئ‎ ٤0 


© = angular direction of upwelling beam (68,9), 

p(Q,02,) = scattering phase function and 

mF, = incoming solar radiative flux. 
The first term on the right hand side describes attenuation of the beam by absorption 
or scattering; the second accounts for intensification by multiple scattering into the 
beam in the following sense. L(t,Q2’) is the diffuse radiation from all directions; 
p(92,02’) accounts for the part of L which has been scattered into the viewing angle, 
2820 ۵/4۲ accounts for the absorption of L scattered into the viewing angle. The third 
term is the intensification by single scattering; that is, the solar irradiance which is 


directly scattered toward the sensor. 
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Fig. 2.4 Relation of scattering, zenith, and azimuth angles, from 
Liou (1980). 
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The multiple scattering term is difficult to evaluate because the intensity field 
must be integrated over all directions. It 1s reasonable to make a single scattering 
approximation if the optical depth is small and the mean free path between aerosol 
particles is large. Global values of optical depth range from 3.0 for a boundary layer 
with heavy haze to 0.003 for the troposphere (World Climate Research Programme, 
1986). Typically, the aerosol optical depth over the ocean is less than 0.2. Therefore it 
is assumed that aerosol particles are randomly distributed in air, the number of 
particles is small and their separation is large, so that single scattering of radiation 
dominates. Also, the upward intensity at the bottom of the atmosphere 1s assumed to 
be zero (L(tj;M,$)=0). The RTE for an atmosphere with total optical depth of t, is 
then: 


L(0;4,@) = (2.8) 


(eon Fo/4(u + no) p(O)L1-exp(- t, CI/B 9 1/09))]), 


where © is the single scattering angle. For an optically thin atmosphere 
t(l/u+1/M )< <1; for example, given t=0.04 and 0,0,= 10° (4,4, > 0.985), 
Tı(l/M ۱/۳۵(< 0۰08, Then, the RTE reduces to: 


L(0:1,9) ^ (© 5F 9/48) p(O)ty, (2.9) 


which says that the reflected intensity from single scattering in an optically thin 
atmosphere is a nearly linear function of tj. Since @,~1, the slope depends on the 
scattering phase function and the satellite viewing geometry. Durkee (1984) calculated 
correlation coefficients for various combinations of terms in equation 2.9 (see Table 1). 
The largest improvement in the coefficient, indicating increased precision, resulted from 
the addition of u, indicating the satellite viewing geometry is the most important 
information for accurately describing the relationship between L and t. 

Table 2 shows the variation of scattering coefficient (s) and single 
scattering albedo (@,) with wavelength for the marine, rural and urban aerosol models 
of Shettle and Fenn (1979). It is to be recalled that Fig. 2.3 shows the variation of 


phase function p(O) with wavelength. It is clear that p and ©, (and, therefore, t) 


scat 
vary more with wavelength for rural and urban aerosols than for marine aerosols. 
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TABLE ! 
CORRELATION ON COMPARISON OF 


FITE 
AND AEROSOL DEPENDENT TERMS IN EQN. 2.9 


Terms compared Correlation coefficient 
L vs t 0.11 
4uL/F, vs t 0.56 
4uL/F, vs p(O)t 0.67 


DE) 


There is little variation in p(O). The result is that €), and © are larger at red-visible 
than at near-infrared wavelengths for rural and urban aerosols; but they are relatively 
constant for marine aerosols. 


Using equation 2.9 for two wavelengths, a ratio can be formed: 


The right hand side will be larger for continental aerosol types (rural and urban) than 

for marine aerosols and, therefore, ratios of radiance at those wavelengths will be also. 
Eqn. 2.10 describes the ratio for aerosol radiance. The total diffuse radiance, 

L, 1s the sum of radiance due to Rayleigh scattering, Lp, and radiance due to aerosol 


scattering, Ly. Therefore the ratio in Eqn. 2.10 can be expressed as: 


Lred/ENIR = (LA + LR)red/(La + LR)NIR: (2.11) 


The relative importance of Rayleigh scattering was estimated by calculating 
Rayleigh radiance values at 0.63um and 0.86um. Climatological April atmospheric 
profiles at 5O°N, the equator, and 50°S were used to estimate variation in Rayleigh 
optical depth (McCartney, 1976). Rayleigh radiance was then calculated for various 
sun-earth-satellite geometries. Table 3 shows the values obtained. Notice how small 
the optical depth variation is with latitude and how little the ratios vary with latitude 
or geometry. 

he ratio of Rayleigh optical depths indicates that if few aerosols are present 
the ratio values should be 3. Since the ratio varies so little in a purely Rayleigh 
atmosphere, the magnitude of (2.11) depends on the spectral variation of La. To 
calculate a correct ratio of aerosol radiance it is therefore necessary to estimate the 
Rayleigh contribution, Lg. In this study the Rayleigh radiance was not explicitly 
removed from the satellite-measured radiance. However, regions with radiance much 
above the Rayleigh radiance (see Table 3) are dominated by aerosol scattering. 
Therefore, where L is significantly larger than Lp the ratio is presumed to be a good 


indicator of aerosol type. 
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TABLE 2 


SCATT 0(۳ ICIENT AND SINGL 


FICIE E SCATTERING ALBEDO 
ROSOL TYPE AND WAVEL Jt 


ENGTH 


50% RI 80% RH 
À Rural Urban Marine 
0.694 0.106 0.091 0.207 
9 , 06 0.050 0.049 0.191 
0.694 0.934 0.635 0.994 

سنا 


1.06 0.865 0.567 0.992 
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C. ALBEDO AND RADIANCE 


Albedo is defined as the ratio of outgoing radiance to incoming radiance; i.e., 
EL (2.12) 


and is usually expressed as a percentage. The ratio of albedos at two different 
wavelengths is directly proportional to the ratio of outgoing radiance and, for the 
optically thin plane parallel atmosphere, to the optical depth NOAA AVHRR 
radiometer brightness counts can be directly transformed into albedos using linear 
relationships derived in pre-flight calibrations (Kidwell, 1984). Therefore, the results 
presented here are in terms of an albedo derived from the NOAA-7 AVHRR 


measurements. 
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TABLE 3 
RAYLEIGH OPTICAL DEPTII, RADIANCE AND RATIOS 


SON 


1.9903 
oid 
2.4783 
1.3150 


Equator 


1.8836 
7770 
2.3475 
1.2439 


SOS 
1.9701 
1.8586 
2.4534 
1.3015 
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III. DATA ANALYSIS 


A. SATELLITE AND DATA SET DESCRIPTION 

The data set used in this thesis was obtained using NOAA-7, a sun-synchronous 
polar orbiting satellite at a nominal altitude of 833 km. It crosses the equator at 
approximately 0230 (ascending; northbound crossing) and 1430 (descending; 
southbound crossing) Local Standard Time. Back scattered radiance is detected with 
the Advanced Very High Resolution Radiometer (AVHRR), a cross-track scanning 
system, in five bands of the solar spectrum called channels. Channels 1 (0.65um), 2 
(0.85um), and 4 (11.0um) were used in this study. The bandwidths of these three 
channels are given in Table 4. The data have a satellite subpoint resolution of 1.1 km. 

Although there is a water vapor absorption band at 0.94j1m and, therefore, in the 
channel 2 band, its effect on the ratio is slight (Durkee, 1984). A decrease in water 
vapor absorption, as would occur over land, would increase channel 2 albedo, and 
therefore decrease the ratio. An increase in water vapor absorption, indicative of a 
more marine environment, would decrease channel 2 albedo, and thus increase the 
ratio. But in a marine environment the increased relative humidity would increase the 
the number of large particles in the size distribution. The channel 2 albedo would then 
increase more than the channel | albedo, thus decreasing the ratio. In this study it is 
assumed that these opposing effects cancel. 

A reduced resolution subset of the AVHRR data was used to do the calculations 
on a VAX11/780 computer. An onboard satellite processor samples four of every five 
picture elements (called pixels) along the scan line to compute an average value. The 
data are further sampled by using every third scan line. This leads to an effective 
subpoint resolution of 4 km by 3.3 km. The particular data used here were copied 
from a larger set at the National Center for Atmospheric Research (NCAR), Boulder, 
Colorado. Data were recorded for six consecutive days, followed by a ten-day lapse, 
and then a second six-day period. The dates were 5-10 and 20-25 April 1982. The data 
set was restricted to the Pacific Ocean (50° N-50°S, 110°E-70°W). Since each pass was 
nominally thirty degrees wide it took seven passes to cover the Pacific Ocean. There 


was some overlap of successive passes. 
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TABLE 4 
AVHRR BANDWIDTHS FOR CHANNELS 1,2 AND 4. 


Channel Bandwidth({tm) 
1 0.58-0.68 
2 0.725-1.10 
4 10.5-11.3 
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The AVHRR was pre-flight calibrated such that a linear relationship between 
brightness counts and percent albedo for both channels 1 and 2 was known. Slope and 
intercept values are part of the source data stream from the AVHRR; thus, it was a 
simple procedure to convert channel | and 2 brightness counts to albedo percentages 


and to form their ratio. 


B. SUNGLINT AND CLOUD DETECTION 

The data were checked in various ways to ascertain that the data would be 
representative of aerosol concentration only, and not contaminated by other effects. A 
particularly spectacular case where a haze layer was visible in channel 1 will be used to 
illustrate these checks. Fig. 3.1 shows the channel 1 image. The NCAR data set had 
more thar one satellite pass on each file. Notice that two passes are appended and the 
image changes three-fourths of the way through. The first part is the mid-Pacific 
Ocean, approximately 175°W from 35°N to 50°S. The last part is the west coast of 
North America from 50°N to 30°N. Sunglint is a serious contaminant of channel | 
and channel 2, but is easily detected. The sunglint is a function of viewing geometry 
and sea state. Since the daylight orbits were near 1430 LST the geometry dictated that 
the glint would always be on the western edge of the image. Thus the sea state 
determined the degree of broadness of the glint region. After visual examination of 
many images, a conservative approach required that the entire western half of the pass 
between 40°N and 10°S be deleted. Fig. 3.2 shows an image with the sunglint region 
deleted. 

The ratio of channel 1 to channel 2 albedo proved useful in two ways. 
Vegetation has a high near-infrared reflectance so upwelling radiance from land should 
be greater at near-infrared than at visible wavelengths. Therefore ratio values over 
land should be smaller than over water, and were in fact less than 1.0. Fig. 3.3 shows 
the west coast of the U.S. deleted, at least where there is no cloudiness. Vancouver 
Island, British Columbia is labelled for reference. This ratio was also one of several 
checks used to identify small cumulus clouds and is discussed further below. 

Three checks were used to screen out clouds. First, channel 2 images were 
visually inspected to determine a threshold value above which it was certain there were 
clouds present. This check was good at eliminating thick (therefore highly reflective) 
low, middle and upper-level clouds, but it did not eliminate thin cirrus or small fair 


weather cumulus. Fig. 3.4 has pixels above the chosen threshold deleted. 
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Fig. 3.1 NOAA-7 unenhanced AVHRR channel | image from 
6 April 1982 covering the central Pacific (~ 175° 
from 50°N to 50°S, and the west coast of North America from 50°N to 30°N. 
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bız Rae as Fig. 3.1 with sunglint region deleted 
40°N to 10°S, western half of image. 
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Fig. 3.3 Same as Fig. 3.1 with sunglint, land and cloud 
pixels deleted. 
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Fig. 3.4 Same as AU 3.1 with sunglint and clouds deleted 
based on channel 2 brightness value threshold test. 
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Second, channel 4 temperatures were used to detect thin cirrus, which is much 
colder than the ocean surface, but is often undetectable in channels 1 and 2. Again, a 
conservative approach was taken. Any channel 4 pixel value greater than 250 was 
eliminated due to possible contamination by cirrus. 

The third check involved situations where a thick aerosol layer resulted in 
channel | and 2 brightness counts at or slightly above brightness counts representative 
of small fair weather cumulus clouds. It was difficult to pick out the cloudy pixels 
since their temperature was close to the ocean temperature. Thus these pixels were 
undetectable in channel 4. The ratio of channel 1 to channel 2 albedo turned out to be 
a good tool in these cases. Small cumulus in clear air had ratios near 1.2 whereas 
aerosol was at or above 1.4. Ratios in clear air over open ocean were around 1.8. 
Therefore 1.3 was chosen as the upper limit for the ratio value, below which it was 
assumed that small clouds were contaminating the scene. Table 5 shows representative 
ratio values for various conditions. These two tests (channel 4 and the ratio values) 
were applied simultaneously so that cloudless cold water regions in the far North and 
South Pacific Ocean would not be excluded by the channel 4 test alone and small 
cumulus, which were missed by channels 2 and 4 alone, would be excluded. In Fig. 
3.5, it is to be noted that the clouds visible in Fig. 3.1 near 50°S are removed and the 
water areas near SO°N are retained. It was very difficult to eliminate the smallest 
cumulus elements based on these tests, especially when they laid under thick haze. 

It is obvious from the preceding figures that each test by itself takes out most 
unwanted pixels. There is much overlap in applying the tests. Yet when applied 
together it is estimated that more than 95% of the unwanted pixels are deleted from 
the image, based on a subjective visual examination of images. The resulting image for 
our case study is Fig. 3.6. The haze layer now stands out clearly. A blow-up of this 


region is shown in Fig. 3.7. 


C. STATISTICS 

The statistics were computed by averaging in one degree boxes over the region. 
This means the area was split into a 181 X 101 array. In addition to mean albedos for 
channels 1 and 2 and their mean ratio, standard deviations of the values were 
calculated. A time average over the first six days (5-10 April), the second six days 


(20-25 April) and the entire data set was done, and the results are presented separately. 
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E CHANNEL I/CHANNEL 2 
IOUS ATMOSPHERIC CONDITIONS 


Ratio value 


0.96 


1.18 


1.20 


1.50 


1.59 


1.78 


TABLE 5 


> 
= 
ہے پر 


Region 
cloud-free land 
cirrus in a haze-free ocean 
cumulus in a haze-free ocean 
small cumulus under haze 
cloud-free haze in occan 


cloud-free, haze-free occan 
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Fig. 3.5 Same as Fig. 3.1 with sunglint and clouds deleted 
based on channel 4 brightness value threshold and 
ratio value threshold tests. 
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Fig. 3.6 Same as Fig. 3.1 with all contaminant tests turned on: 
this image was included in climatological statistics. 
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E 3.7 Enlargement of a portion of Fig. 3.6: 
right strip across center of image 1s haze. 
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The number of “good” pixels (not cloudy and not in the sunglint over open 
ocean) was retained for each time period. Therefore in each one degree box it was 
known how many pixels were used in the statistics. Then, those areas of the North 
and South Pacific Ocean which had better data coverage, and therefore where the 
statistics were more reliable, could be estimated. This also gave an indication of where 
persistent cloudiness occurred. 

It can be seen from Fig. 3.8 that the sunglint check between 40°N and 10°S, seen 
as generally northwest-southeast oriented strips, deleted a large number of pixels 
compared to the southern hemisphere. It can also be seen that adjacent overlapping 
satellite passes did not fill in these strips evenly. This is an artifact of the small data 
set used and would presumably disappear with a sufficiently large sample. Fig. 3.9 
illustrates the data density for the second period. Although there is much less data 
than for the first six-day period, the same pattern is observed. It can be seen from the 
full data set coverage pattern, Fig. 3.10, that there is a dramatic increase in data south 
of 10°S, the sunglint strips are discernable but more evenly filled in than for the two 
six-day periods, and there is a paucity of data in northern and southern latitudes. 

The major feature in Fig. 3.11, the channel I albedo for 5-10 April, is the 
maximum in the northern hemisphere and the minimum in the southern hemisphere. 
There is a minimum in the Sea of Japan and along the west coast of North America, 
and a maximum on both sides of Central America; one stretching north and then east 
across the Gulf of Mexico, the other to the southeast. The pattern in Fig. 3.12, the 
channel 1 albedo for 20-25 April, is similar, with the maximums near Central America 
expanded. Fig. 3.13 illustrates the resulting pattern for the full data set. It can be seen 
that there is a maximum in channel | albedo in the northern hemisphere and a 
minimum in the southern hemisphere, with distinct regional minima in the Sea of 
Japan, and along the west coast of North America, and regional maxima on both sides 
of Central America. Recall from the discussion of Table 3 that there is a critical value 
for the albedo, above which it is presumed aerosol scattering is dominating Rayleigh 
scattering; that value is approximately 2.0. Therefore, in Figs. 3.11, 3.12 and 3.13, 
where the values of albedo are significantly above 2.0, it is presumed that aerosol type 
can be inferred from the ratios of channel 1 to channel 2 albedo. Note that this 1s true 
for nearly all of the area. 

The El Chichön volcano (17.3°N, 95.2°W) erupted several times between 28 
March and 6 April, spewing debris into the troposphere and stratosphere. Robock and 
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Fig. 3.8 Number Pee used in calculating statistics for 
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Fig. 3.10 Number of D used in calculating statistics for 
5-10 and 20-25 نک‎ 1982. Areas with greater than 
pixels are outlined. 
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1 Mean channel 1 albedo for 5-10 April 1982. 
Albedos greater than 6.0% are outlined. 
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Mean channel | albedo for 20-25 April 1982. 


Albedos greater than 6.0% are outlined. 
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Fig. 3.13 Mean channel I albedo for 5-10 and 20-25 April 1982. 
Albedos greater than 6.0% are outlined. 
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Matson (1983) mapped the circumglobal transport of the stratospheric cloud in the 
months after the eruptions, based on GOES and NOAA-7 visible and infrared images. 
Krueger (1983) estimated the positions of both stratospheric and tropospheric clouds in 
the days immediately following the eruptions based on Nimbus-7 total ozone mapping 
spectrometer data. He indicated the tropospheric debris cloud moved west while the 
Stratospheric cloud moved east, based on observed winds from rawinsonde ascents. 
Sea level pressure and 500 mb analyses for 12 GMT 5 April and 12 GMT 20 April 
1982 were referred to in an attempt to verify that the albedo patterns seen in Figs. 3.11, 
pl? and 3.15 could) be ascribed to the 2] Chichon eruptions. Figs. 3.14 and 3.15 
support the supposition that the El Chichon cloud is detectable in the NOAA-7 
channel | albedo patterns of Figs. 3.11, 3.12 and 3.13. 

Fig. 3.16 illustrates the channel 1 albedo standard deviation pattern for 5-10 
April. It can be seen that maximum standard deviation occurs in the tropics. Fig. 3.17 
illustrates the channel | albedo standard deviation pattern for 20-25 April. The pattern 
again shows highest standard deviation in the tropics. It can be seen from Fig. 3.18, 
the full data set channel 1 albedo standard deviation, that values are highest near the 
equator, and lowest near the extreme northern and southern latitudes, and along the 
west coast of North America. The reasoning for this pattern is not understood. 

The albedo pattern in Fig. 3.19 for channel 2 during 5-10 April shows all of the 
same features as in channel | albedo (compare Fig. 3.11), except the albedo magnitude 
is smaller for channel 2. The northern hemisphere is a maximum, while the southern 
hemisphere 1s a minimum. There are minima in the Sea of Japan and along the North 
American west coast, and maxima on either side of Central America. Fig. 3.20 
illustrates the channel 2 albedo pattern for 20-25 April. The pattern is again similar to 
channel 1 albedos for the same period (compare Fig. 3.12). The channel 2 albedo 
pattern for the full data set, illustrated in Fig. 3.21, is remarkably similar to the channel 
1 albedo pattern (compare Fig. 3.13). Maximum albedos occur in the northern 
hemisphere; minimum albedos occur in the southern hemisphere. The regional maxima 
and minima seen in the two six-day periods are evident here. Channel 2 albedos are 
smaller in magnitude than channel 1 albedos. 

Channel 2 albedo standard deviation for 5-10 April is shown in Fig. 3.22. 
Maximum values occur in equatorial regions; minimum values occur in the extreme 
latitudes and along the west coasts of North and Central America. Fig. 3.23 illustrates 


the channel 2 albedo standard deviation for 20-25 April. It shows a maximum in the 
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Fig. 3.14 Sea level pressure analysis for 12 GMT 05 April 1982. 
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Fig. 3.15 500 mb analysis for 12 GMT 05 April 1982. 
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Fig. 3.16 Channel 1 albedo standard deviation for 5-10 April 1982. 
Standard deviations greater than 1.5% are outlined. 
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Standard deviations greater than 1.5% are outlined. 
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Fig. 3.18 Channel 1 albedo standard deviation for 5-10 and 20-25 
April 1982. Standard deviations greater than 1.5% are outlined. 
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Fig. 3.19 Mean channel 2 albedo for 5-10 April 1982. Albedos less 
than 2.0% and greater than 4.0% are outlined. 
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Albedos 


less than 2.0% and greater than 4.0% are outlined. 
54 


Fig. 3.20 Mean channel 2 albedo for 20-25 April 1982. 
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Fig. 3.21 Mean channel 2 albedo for 5-10 and 20-25 April 1982. 
Albedos less than 2.0% and greater than 4.0% are outlined. 
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tropics, although not as extensive as for 5-10 April. It can be seen from the channel 2 
albedo standard deviation pattern for the full data set, shown in Fig. 3.24, that 
maximum values occur near the equator, whereas minimum values occur at the 
extreme latitudes and along the west coasts of North and Central America. 

Recall that it was difficult to remove all small cumulus clouds from the data base. 
This does not present a serious problem because in the marine environment, clouds 
cause a bias toward larger particle size; in other words, the aerosol would appear more 
marine in character. Ratio values would be decreased by the presence of small 
cumulus clouds. Having acknowledged this bias, the channel 1 to channel 2 albedo 
ratio pattern for 5-10 April can be seen in Fig. 3.25. There is a maximum in the 
tropics, with regional maxima in the East China Sea, the Sea of Japan and east of 
Japan, off the north, west, and south coasts of Australia, around New Zealand, and off 
the west coasts of Mexico and South America. Minimum values occur at the extreme 
latitudes, with regional minima in the southwestern Gulf of Mexico, off the east coast 
of Australia, and southeast of Japan. Fig. 3.26 illustrates the ratio pattern for 20-25 
April. Maximum values occur in the tropics, with regional maxima in the Sea of 
Japan, off the northwest, south and southeast coasts of Australia, north of New 
Zealand, and off the west coasts of North, Central and South America. Minimum 
values are found at the extreme latitudes, and in the Gulf of Mexico, east of Australia 
and southeast of Japan. Ratio values for the full data set are illustrated in Fig. 3.27. 
Maximum values extend from the west coasts of North, Central and South America 
across the equatorial Pacific Ocean. Maxima can also be seen around and northeast of 
New Zealand, in the East China Sea and Sea of Japan, and off the north, west and 
south coasts of Australia. Minimum values occur in extreme northern and southern 
latitudes, the southwest Gulf of Mexico, southeast of Japan, and off the east coast of 
Australia. 

Recall that, if an estimate for Rayleigh scattering can be made, high ratio values 
indicate the presence of a continental type of aerosol and low ratio values indicate a 
marine aerosol distribution. It can be seen from Fig. 3.27 that China, North America 
and Central America are apparent sources of continental aerosols in the northern 
hemisphere, and Australia, New Zealand and South America are apparent sources of 
continental aerosols in the southern hemisphere. Sea level pressure and 500 mb 
analyses were referred to in an effort to support this interpretation of the data. Surface 


winds over northern Australia were light and variable in the early days of April, while 
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Fig. 3.22 Channel 2 albedo standard deviation for 5-10 cub 1982. 
Standard deviations greater than 1.5% are outlined. 
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3.23 Channel 2 albedo standard deviation for 20-25 April 
982. Standard deviations greater than 1.5% are outlined. 
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Fig. 3.24 Channel 2 albedo standard deviation for 5-10, 20-25 
April 1982. Standard deviations greater than 1.5% are outlined. 
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Fig. 3.25 Mean channel 1/channel 2 ratios for 5-10 April 1982. 
Ratios greater than 1.6 are outlined. 
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for 20-25 April 1982. 
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nel 1/channel 2 r 
1982. Ratios greater than 1.6 are outlined. 
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two tropical storms formed to the northeast. Later in the month easterly trade winds 
re-established themselves over the region. The west and south coasts of Australia were 
influenced by a high pressure cell that moved from the Bight to Tasmania during the 
month. This kept the winds off-shore on those coasts. Eastern Australia, as well as 
New Zealand, were affected by the two tropical storms. Both storms moved 
southeastward, the first becoming an intense extra-tropical storm that passed just to 
the north of New Zealand, the second slowly dying off the east coast of Australia. 
These systems kept the winds on-shore in eastern Australia and brought southerly 
winds of 40 kt to New Zealand on 9 and 10 April. The remainder of the month 
brought southeasterly winds of 15-25 kt to New Zealand. At 500 mb, winds over 
Australia were generally westerly and increased through the month. Over New Zealand 
500 mb winds were west-southwesterly. 

In the China-Japan area, a strong high pressure cell at sea level in early April 
moved south and weakened, allowing several lows to develop near southern Japan and 
move rapidly northeast. These storms brought strong west-northwesterly winds out of 
China and across the East China Sea and Sea of Japan while the trade winds south of 
Japan remained strong out of the east. 

The apparent maximum of continental aerosol throughout the tropical Pacific 
© is the most unexpected feature in the ratio pattern. The convergence of 
northern hemisphere northeasterly trade winds and southern hemisphere southeasterly 
trade winds, carrying continental aerosols from the Americas, is a possible explanation 
of the equatorial Pacific Ocean maximum. Time sequences and trajectory analyses 
must be studied to fully understand this phenomenon. 

The Gulf of Mexico minimum in ratio values in the region of El Chichón is at 
first surprising. However, Dutton and DeLuisi (1983) have shown that the spectral 
response ~f the aerosol optical depth for the El Chichón cloud was very flat (see Fig. 
3.28). A ratio of channel 1 to channel 2 albedos from an aerosol with optical depth as 
in Fig. 3.28 would result in low ratios. 

Fig. 3.29 illustrates the ratio standard deviation pattern for 5-10 April. 
Maximum values occur in the tropics and East China Sea; minimum values in the 
southern hemisphere, southwest of North America and in the Gulf of Mexico. The 
20-25 April ratio standard deviation pattern is shown in Fig. 3.30. The pattern is 
similar to that in Fig. 3.29, but maximum values are not as high in the tropics. [t can 


be seen from the full data set ratio standard deviation pattern of Fig. 3.31 that 
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maximum values occur in the equatorial Pacific Ocean and East China Sea. Minimum 
values occur in the southern hemisphere, northern hemisphere except for the northwest 


Pacific Ocean, and southwestern Gulf of Mexico. 
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Fig. 3.28 The El Chichon co وت‎ of aerosol optical depth as a 
function of wavelength. Curves are labelled with station 
or latitude and ee (from Dutton and DeLuisi, 1983). 
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Fig. 3.29 Ratio standard deviation for 5-10 و‎ 1982. 
Standard deviations greater than 0.2 are outlined. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A NOAA-7 AVHRR data set for 5-10 and 20-25 April 1982 was used to establish 


a physical basis for an aerosol climatology in the North and South Pacific Ocean. 


Channel 1 (0.63u0m), channel 2 (0.864m), and channel 4 (119m) data were used to 


eliminate pixels that were over land or cloudy and to compute various statistics. A 


summary of conclusions follows: 


In addition to using threshold values in channel 2 and 4 albedos to eliminate 
clouds, a new technique based on the ratio of channel I/channel 2 albedo was 
developed to eliminate small fair weather cumulus that would otherwise have 
been included in the data set. These cloud elements are difficult to differentiate 
because their brightness counts 1n channel 2 data are cqual or greater than that 
for heavy haze events. Their temperatures from channel 4 data are nearly equal 
to the Sea-surface temperature. However, ratio values were less thån 1.3, 
whercas haze events had ratios of 1.4 or greater. Therefore application of this 
method includes cases of heavy haze where brightness counts would be greater 
than the cloud detection threshold. 


High ratio values were obtained over the equatorial Pacific Ocean, along the 
west coasts of North, Central and South America, 1n the Sea of Japan and East 
China Sea, off the north, west and south coasts of Australia, and around New 
Zealand. Low ratio values were obtained in the extreme northern and southern 
latitudes, southeast of Japan, the southwestern Gulf of Mexico, and off the east 
coast of Australia. Theory leads to the conclusion that high ratio values 
indicate areas of continental aerosol type, while low ratios indicate areas of 
marine aerosol type. This, in, turn, suggests that Asia and North and Central 
America are sources of continental aérosol in the northern pe and 
Australia, New Zealand and South America are sources of continental aerosol 
in the southern hemisphere. 


Channel 1 and 2 albedos can be used to identify major aerosol-producing 
events, such as volcanic eruptions. A pattern in the albedos was linked to the 

Chichon eruptions of March-April 1982, based on previous analyses of the 
debris cloud drift and winds from surface and upper-air data at the time. 


Several improvements and further studies remain to be done. The cloud 


l. 


یا 


discrimination technique needs to be refined to ensure a complete elimination of clouds 


while keeping all cloud-free regions in the data set. The sun-earth-satellite geometry 


can be taken into account, correcting the cloud detection threshold for the latitudinal! 


variation in cloud brightness. The effects of water vapor absorption on AVHRR 


channel I and 2 albedos, and, therefore, on ratio values, should be explicitly studied. 


The contribution of Rayleigh scattering to the total atmospheric scattering should be 


accounted for in the determination of ratio values. Anomalously high ratio values were 


found in several small areas of the study basin, both over land and at sea. Values from 


5 to as high as 20 were seen. These should be investigated for physical explanations. 
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After a larger data base is available, time series and trajectory analyses can be 
done to investigate aerosol dispersion and seasonal variation. The mid-Pacific Ocean 
equatorial maximum in ratios also needs to be more closely investigated. Finally, 
further determination of sources of continental aerosols should be made. 

This study represents the first step toward a satellite-derived aerosol climatology. 
Once fully developed, this method for establishing climatologies of aerosol amount and 
particularly aerosol type will be a much needed addition to studies of atmospheric 


radiation, the general circulation and long-term climate variations. 
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